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FOREWORD 


The  work  in  this  report  is  prepared  for  the  Army  Materials  Technology  Lab¬ 
oratory  (MTL)  under  Contract  No.  DAAG46-85-K-001 1  with  the  Institute  of  Fracture 
and  Solid  Mechanics  at  Lehigh  University,  Bethlehem,  Pennsylvania.  The  authors 
are  particularly  grateful  to  the  Project  Manager,  Mr.  J.  F.  Dignam  and  Technical 
Monitor,  Dr.  S.-C.  Chou  from  MTL.  Their  constant  interest  and  encouragement  have 
made  the  progress  in  this  work  possible. 

An  improved  version  of  the  computer  code  is  made  for  predicting  the  damage 
of  projectile/target  systems  where  failure  takes  place  predominantly  by  perma¬ 
nent  deformation,  fragmentation  and/or  fracture.  Instead  of  considering  the 
failure  of  individual  elements,  a  more  detailed  account  of  the  damage  pattern  can 
be  obtained  from  the  distribution  of  the  surface  and  volume  energy  density.  In 
order  to  demonstrate  the  validity  of  the  energy  density  theory,  it  is  only  logi¬ 
cal  to  first  consider  an  example  in  which  the  projectile  velocity  is  sufficiently 
low  such  that  material  phase  transformation  does  not  occur.  The  method  of  ap¬ 
proach,  however,  applies  equally  well  to  situations  where  the  solid  can  transform 
to  liquid  and/or  gas  as  the  projectile  velocity  is  increased  into  the  hyperveloc¬ 
ity  range.  Unlike  the  conventional  approaches,  the  energy  density  theory  requires 
only  a  knowledge  of  the  initial  uniaxial  response  of  the  projectile/target  system. 
Change  in  the  local  strain  rates  are  determined  analytically  for  each  time  incre¬ 
ment.  The  proposed  methodology  has  far-reaching  significance  as  it  can  provide 
useful  information  for  modeling  subscaling  tests. 
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INTRODUCTION 


Damage  prediction  of  projectile/target  systems  has  been  the  subject  of  many 
past  discussions  [1-4].  The  majority  of  the  works  make  use  of  the  classical 
theories  of  continuum  mechanics  by  assuming  that  the  constitutive  relations  of 
the  materials  are  specified  or  known  as  an  a'  priori.  These  theories  are  general¬ 
ly  referred  to  as  elasticity,  plasticity,  viscoplasticity,  etc.,  and  are  not  suf¬ 
ficient  for  describing  the  projectile/ target  damage  process  unless  additional 
failure  criteria  are  invoked.  Such  a  loosely  adopted  procedure  leaves  room  for 
arbitrariness  in  terms  of  empirical  parameters  and  lacks  the  consistency  and 
uniqueness  required  in  theoretic  prediction.  Stress  and  failure  analysis  cannot 
be  separated.  They  should  be  inherently  coupled  in  the  theoretical  formulation. 

An  earlier  attempt  was  made  in  [5]  to  analyze  the  damage  of  a  project! le/tar- 
get  system  by  application  of  the  energy  density  theory  [6,7].  Data  banks  were 
constructed  to  account  for  the  change  in  material  behavior  of  the  projectile  and 
target  as  damage  occurred  progressively.  Changes  in  local  strain  rates  and 
strain  rate  history  are  included  such  that  the  dynamic  stress  and  stress  behavior 
of  each  material  element  is  different  and  derived  independently.  For  the  first 
time,  stress  and  failure  analysis  were  treated  simultaneously.  Damage  by  perma¬ 
nent  deformation  and/or  fragmentation  was  modelled  by  application  of  triangular 
finite  elements.  Element  failure  is  assumed  to  occur  when  the  corresponding  sur¬ 
face  and/or  volume  energy  density  reach  their  critical  values.  These  elements 
are  thus  removed  from  the  analysis  and  the  grid  pattern  is  modified  accordingly. 
This  process  is  repeated  for  each  time  increment  until  the  target  is  completely 
fractured.  The  shortcoming  of  the  element  failure  model  became  obvious  as  the 
solution  accuracy  depended  exclusively  on  the  size  of  the  finite  element  mesh. 


Increasing  accuracy  can  only  be  achieved  at  the  expense  of  extensive  labor  and 
computer  time.  Nevertheless,  the  qualitative  feature  of  target  failure  by  plug¬ 
ging  is  predicted.  What  has  been  accomplished  is  a  preliminary  confirmation  on 
the  validity  of  the  methodology. 

An  alternative  scheme  for  determining  the  progressive  damage  of  projectile/ 
target  systems  is  developed  in  this  investigation.  Instead  of  considering  the 
failure  of  individual  elements,  contours  of  surface  and  volume  energy  density 
are  obtained.  Damage  patterns  are  then  mapped  in  accordance  with  the  intensity 
distribution  of  the  energy  densities.  This  may  involve  the  creation  of  free 
surface  areas,  fragmentation  over  a  region  or  the  fracture  of  solid  segments. 

A  much  more  realistic  prediction  of  projectile/target  failure  behavior  is  ob¬ 
tained.  Even  though  results  are  displayed  only  on  the  mechanical  fragmentation 
and  plugging  of  the  impacted  target,  change  in  local  temperature  and  the  forma¬ 
tion  of  adiabatic  shear  band  [8-10]  are  also  contained  in  the  present  work.  Ma¬ 
terial  transformation  that  occurs  within  this  narrow  band  of  material  adjacent 
to  the  plug  can  be  verified  by  computing  for  the  latent  heat.  Other  thermal 
properties  of  the  projectile/target  material  may  be  made  available  to  check  with 
available  experimental  data.  Study  on  changes  in  microstructure  is  beyond  the 
scope  of  this  investigation. 

PROJECTILE/TARGET  DAMAGE  CONSIDERATION 


Damage  is  a  process  by  which  material  elements  experience  permanent  change 


and  no  longer  return  to  their  original  states.  It  changes  with  time  and  continues 
to  occur  if  the  surface  and  volume  energy  density  are  not  evenly  proportioned. 

As  a  result,  the  rate  change  of  volume  and  surface  area  of  each  element,  i.e., 


dV/dA  depends  on  the  rate  of  energy  input.  The  energy  quantities  identified 


with  dV  and  dA  are,  respectively,  the  volume  energy  density  dW/dV  and  surface 


energy  density  dW/dA.  Their  interdependence  is  well-known  in  the  homogeneous 


nucleation  theory  of  Gibbs.  The  inclusion  of  the  length  parameter  dV/dA  in  the 


energy  density  theory  represents  a  fundamental  departure  from  all  existing  ap¬ 


proaches, 


Vamage  Plana.  Each  unit  volume  of  material  in  the  projectile/target  system 


is  considered  to  be  damaged  under  deformation  or  strain,  say  If  the  rec¬ 


tangular  Cartesian  coordinates  (x,y,z)  are  used  to  denote  the  original  coordi¬ 


nates,  then  the  deformed  coordinates  locate  the  so-called  "plane  of  homo¬ 


geneity"  or  damage  plane.  An  unique  damage  or  strain  state  is  said  to  exist  for 


each  time  increment  of  loading.  It  is  only  on  this  plane  that  the  volume  en¬ 


ergy  density  can  be  completely  determined  by 


The  quantity  A  is  a  proportionality  parameter  that  relates  the  slope  of  the 


equivalent  uniaxial  stress  and  strain  relation  on  the  damage  plane  to  the  rate 


of  change  of  volume  with  the  surface  area,  i.e.. 


de.  ^^dA^C 


The  surface  and  volume  energy  density  are  assumed  to  be  unconnected  in  the 
classical  theories  of  continuum  mechanics  because  of  the  invokernent  of  dV/dA 
0  in  the  limit.  Such  a  simplification  leads  to  an  inaccurate  account  of  ener¬ 
gy  transmission  across  the  projectile/target  interface. 


The  invariant  concept  so  widely  used  for  deriving  constitutive  relations  in 
continuum  mechanics  no  longer  holds  here  because  the  damage  plane  occupies  a  defi¬ 
nite  orientation  and  is  different  for  each  material  element. 


Hence,  the  surface  energy  density  (dW/dA)  can  be  written  as 


MA^^  dV 


(3) 


with  dW/dV  being  a  scalar.  For  a  homogeneous  and  isotropic  material,  the  dam¬ 
age  plane  can  be  determined  from  the  assumption  that  there  prevails  three  orthog 
onal  directions  (^,n,s)  along  which  the  surface  energy  densities  (dW/dA)^.  (i 
=  C,n,?)  are  equal ,  i .e. , 


I dW \  _  / dW \  _  / dW \ 

~  ^dA^n  ■ 


dA'o 


(4) 


This  provides  an  one-to-one  correspondence  between  ea..h  element  in  the  projec¬ 
tile/target  system  under  a  multi -axial  stress  or  strain  state  with  that  in  the 


uniaxial  data  bank  or  the  value  (dW/dA)^.  Since  equation  (3)  may  be  generalized 


to  read  as 


^dA  i  ”  'Hfl''-' 


^dAM  dV’ 


1  =  C>n,c 


(5) 


it  follows  also  that  the  (dV/dA).,  (dV/dA)  and  (dV/dA),  are  equal  on  the  damage 

.  n  4 

plane. 

Change  VcCaive  a'lth  Aaea.  The  quantity  dV/dA  is  most  unique  as  it  is  re¬ 
lated  to  the  slope  of  the  uniaxial  stress  and  strain  curve  as  indicated  in  equa¬ 
tion  (2).  For  a  two-dimensional  stress  state  in  Figure  1(a),  the  corresponding 
plane  of  homogeneity  or  damage  plane  makes  an  angle  a  with  the  x-  and  y-axis  as 
shown  in  Figure  1(b).  The  stresses  and  strains  are  such  that  the  condition  in 
equation  (4)  is  satisfied.  In  two  dimensions,  the  damage,  say  in  the  z-  or  c- 


I 


direction,  may  be  regarded  as  constant.  This  means  that  (dW/dA)  or  (dW/dA)^ 
remains  the  same.  Alternatively,  it  can  be  stated  that 


(^)  or  (^)  =  c 

^dA^z  °  ^dAV  ^ 


(6) 


where  c  is  a  constant.  Without  going  into  details  and  making  use  of  the  results 
in  [7],  the  strain  component  in  the  z-  or  c-direction  may  be  written  as 


:  or  c,  =  (c-1  ){z^+z  )cos.3 


7) 


It  follows  that  the  (dV/dA)^.  (i  =  c,n)  expressions  are  given  by  [7] 


A 

^dA^ 


A 

ma7 


[(^)s,  +  c  Icosa 


(8) 


£  +£ 


This  is  referred  to  as  the  planar  formulation  of  the  energy  density  theory  [6]. 
The  case  when  c  -«■  1  corresponds  to  the  state  of  plane  strain  in  the  classical 
theory  of  elasticity  with  -  0.  In  general,  c  can  take  any  appropriate  con¬ 
stant  greater  than  unity  depending  on  the  constraint  in  the  :-di recti  on.  The 
determination  of  the  various  quantities  will  be  discussed  in  connection  with  the 
finite  element  procedure. 


Damage  T/iaei/ictdi .  Damage  of  the  projectile/target  system  will  be  analyzed 
by  observing  the  contours  of  dW/dA  on  the  damage  plane  and  dW/dV.  Initially,  the 
surface  energy  density  would  be  more  dominant  as  the  projectile  first  hits  the 
target  where  energy  is  transferred  across  the  contact  or  interface.  Fragmenta- 


-6- 


tion  may  occur  as  the  result  of  (dW/dA)^  exceeding  its  critical  value.  Increase 
in  the  volume  energy  density  takes  place  when  the  disturbance  spreads  over  a 
larger  portion  of  the  target.  The  following  conditions  may  serve  as  a  guide  for 
determining  the  damage  pattern: 


dv  ^ 


dV'c’  ^dA^; 


/CWx 

^dA^c’  ^dA^?  ^dA^c 


(9) 


or 


,dW, 

'dA^ 


MA^c’ 


dW 

dV 


(^)  A  >  A 

MV^c’  ^dA^-:  M.A^c 


:io; 


The  other  possibility  would  be  for  both  (dW/dA)^  and  dW/dV  to  exceed  their  criti¬ 
cal  values  at  the  same  time.  This  would  occur  when  damage  occurs  at  large. 
Equations  (9)  and  (10)  suffice  for  projectile/target  systems  where  failure  is 
dominated  by  permanent  deformation,  fragmentation  and  fracture  including  the 
formation  of  a  plug.  Predicted  damage  patterns  are  simulated  by  altering  the 
connectivity  of  the  triangular  elements  as  displayed  in  Figures  2(a)  and  2(b). 
Regridding  is  then  carried  out  as  further  damage  occurs  within  the  system.  This 
procedure  is  continued  uncil  the  projectile  either  penetrates  through  the  target 
or  arrest  in  the  target. 

Should  the  condition  in  equation  (9)  be  met  locally  over  a  relatively  narrow 
region  such  as  those  shown  between  elements  A  and  B,  C  and  D  and  E,  F  in  Figure 
2(a),  the  enclosed  material  may  be  regarded  as  fragmented  into  small  pieces  and 
their  masses  are  redistributed  onto  the  neighboring  nodes.  The  shaded  areas  are 
then  modelled  as  cavities.  Fragmented  sizes  are  even  smaller  if  the  condition 
(dW/dA),  <  (dW/dA)  prevails  in  equation  (9)  because  (dV/dA)__  would  then  be  much 
smaller  than  (dV/dA)^.  When  the  condition  in  equation  (10)  is  met  by  a  set  of 


constant  (dW/dA).  >  (dW/dA)  contours  that  are  narrowly  spaced  but  encircling 
a  relatively  large  undamaged  area,  then  a  block  of  mass  is  assumed  to  be  dis¬ 
lodged.  lio  damage  is  predicted  in  the  detached  mass  if  dW/dV  <  (dW/dV)^  while 
some  damage  may  be  present  for  the  case  when  dW/dV  -  (dW/dV)^.  The  situation 
just  described  is  shown  in  Figure  2(b)  where  the  element  D  is  dislodged  from 
its  neighboring  elements.  The  element  D  retains  its  stiffness  and  is  trapped 
within  the  system  as  damage  proceeds.  Plugging  would  be  predicted  in  the  same 
way.  The  (dW/dA).  >  (dW/dA)^  lines  would  then  conform  with  plug  shape  in  the 
target.  To  be  emphasized  is  that  (dW/dA).  plays  a  more  dominant  role  in  material 
fracturing  into  block  masses  while  dW/dV  in  material  failure  by  fragmentation. 

METHOD  OF  SOLUTION:  FINITE  ELEMENT 

The  projectile/target  system  is  discretized  by  subdividing  both  of  the  con¬ 
tacting  solids  into  a  finite  number  of  subdomains  or  elements  in  the  form  of  tri¬ 
angles  that  are  interconnected  by  nodes.  A  typical  triangular  element,  say  A  in 
Figures  2(a)  or  2(b),  will  contain  nodes  i,  j  and  k  referred  to  in  a  counterclock 
wise  direction.  The  coordinates  of  each  node  are  located  in  the  xy-domain  such 
that  attention  is  focused  on  analyzing  the  kinematics  of  the  individual  elements. 
A  stiffness  matrix  containing  the  properties  of  both  the  projectile  and  target  is 
obtained  for  relating  the  accelerations  of  each  element  to  the  dynamic  forces  via 
the  nodal  points.  The  incremental  strains  are  then  found  from  the  incremental 
displacements.  As  mentioned  before,  no  a  priori  assumption  is  made  on  the  form 
of  the  constitutive  relations.  What  needs  to  be  known  is  only  the  initial  mate¬ 
rial  behavior  for  the  first  load  increment  which  is  taken  to  be  sufficiently 
small  such  that  all  elements  are  assumed  to  respond  in  the  same  way.  Deviations 
in  the  dynamic  stress  and  strain  behavior  from  element  to  element  are  accounted 
for  thereafter.  There  is  no  restriction  on  the  rate  of  local  strain  which  can  be 


as  high  as  10  or  10  sec  depending  on  velocity  at  which  the  projectile  impacts 
the  target. 

Vodat  .  Referring  to  a  rectangular  Cartesian  coordinate  system 

(x,y),  the  displacement  vector  ^u}  possesses  the  components  u  and  u  and  it 

X  y 

takes  the  form 


In  the  case  of  triangular  elements,  the  displacement  is  a  linear  function  of  the 
space  variables: 


(u)  =  {ci.|}  +  {a2}x  +  {ci3}y 


in  which  {a.|},  {a2}  and  {a^}  are  column  vectors.  It  is  convenient  to  express 
the  displacements  in  terms  of  those  at  nodal  points  {u®}  by  making  use  of  the 
interpolation  function  [N],  i.e.. 


fu)  =  [N]{u®} 


In  equation  (13),  the  nodal  point  displacement  vector  is 


r  e,  r  1  1  J  J  k  k.T 

{u  }  =  (u  ,  u  ,  u  ,  u  ,  u  ,  u  } 

“  .x’  .y’  “x’  ,y’  .X’  “y 


while  [N]  stands  for 

N.,  0  ,  N.,  0  ,  N,  ,  0 
[N]  =  ^  J  " 

0  ,  N.,  0  ,  Nj.,  0  ,  N 


The  components  N^.  are  given  by 


I,  =  jj  (a,  +  b.x  +  c^y) 


with  A  being  the  area  of  the  triangular  element  i,  j  and  k,  i.e., 


1.  X,,  y, 


A  =  2  det  1 ,  X j ,  y. 

"k’  ^k 

The  parameters  a^. ,  b^.  and  c^.  are 


"i  =  ’‘/k  -  Xki'j’  -  l-k’  S  =  Xk  -  Xj 


Once  {u}  or  {u  }  is  known,  the  strains  follow  immediately  from  the  strain/dis¬ 
placement  relations. 

Inciemental  St^icun.  Let  {e}  denote  the  strain  tensor  referred  to  the  com¬ 
ponents  e^,  £y  and  i.e., 


{e}  =s  e 


If  {e}  stands  for  the  strain  at  time  t,  then  the  strain  increment  in  the  time  in¬ 
terval  At  can  be  written  as 


^^f^At  =  ^e^t  ■  ^^^t-At 


A  zero  initial  strain  field  is  taken  while  the  strain  increment  Ae  can  be  ob¬ 
tained  from  Au  as 


■r  -V-";  ■ 


In  terms  of  the  displacement  increments  at  the  nodal  points,  equation  (21)  be¬ 


comes 


{Ae}  =  [B]{Au®} 


(22) 


where 


{AU®}  =  (AU^,  AUy, 


(23) 


The  matrix  [B]  is  defined  by 


[B] 


1_ 

ay  ’  ax 


N.,  0  ,  N^., 
0  ,  N.,  0  , 


0  ,  N^,  0 
N.,  0  , 


(24) 


This  completes  the  preliminaries  on  expressions  associated  with  strain  increments. 


EqLLcvalenX  Uniax^lal  Stteii  Sto-te.  One  of  the  unique  features  of  the  energy 

density  theory  lies  in  the  translation  of  uniaxial  data  to  multi-axial  stress 

or  strain  states  without  loss  in  generality  .  This  is  accomplished  via  the  con- 

*The  plasticity  theory  assumes  the  coincidence  of  the  uniaxial  data  with  the  ef¬ 
fective  stress  and  strain  curve  which  tends  to  ignore  the  effects  of  dilatation 
that  are  by  no  means  negligible  in  the  projectile/target  damage  process. 


cept  of  a  damage  plane  on  which  the  energy  transferred  across  a  unit  surface 
area  is  the  same  in  all  directions.  In  two  dimensions,  such  a  plane  is  defined 


by  [(dV/dA)^,  or  [(dV/dA)^,  which  is  related  to  the  uniaxial  data 
[(dV/dA)^,  e^].  In  view  of  equation  (2),  the  stress  components  7^  and 
tainable  as 


are  ob 


(25) 


Hence,  the  association  of  a.  with  e.  or  7^  with  ■:  may  be  regarded  as  the  equiva 

->  ->  '  •  ^ 

lent  uniaxial  stress  and  strain  relation  because  they  refer  to  an  unique  volume 


energy  density 


dW 

dV 


/  a,d£_  or  r  a  de 
n  n 


(26) 


which  is  the  equivalent  of  equation  (1). 

In  retrospect,  it  can  be  concluded  from  equations  (8)  and  the  condition 

(dV/dA)  =  (dV/dA)  that  e  equals  to  e  on  the  damage  plane.  Equation  (26)  thu 
s  n  s  ^ 

★★ 

follows.  The  equivalence  of  and  however,  is  not  a  necessary  condition 
in  the  energy  density  theory.  It  will  no  longer  hold  when  the  expressions  of 

(dV/dA)^  and  (dV/dA)  in  equations  (8)  are  altered  to  include  nonlinear  and/or 

n 


The  expression  dW/dV  =  /  a^jde^.j  used  in  nonlinear  elasticity  or  plasticity  in 

vokes  the  independence  of  load  sequence  on  stored  energy.  This,  of  course,  is  a 
restriction  that  cannot  be  justified  in  any  theories  involving  irrecoverable  en¬ 
ergy  dissipation. 


This  should  not  be  confused  with  the  plane  of  principal  strain.  Such  a  rela¬ 
tion  occurred  by  chance  as  a  result  of  the  neglect  of  nonlinear  and/or  shear 
stress  effects  in  equations  (8). 


Figure  3.  Schematic  of  actual  stress  and  strain  response  of  a  typical  element. 

in  which  o  is  the  mass  density  of  the  local  element  that  may  vary  from  location 
to  location  in  the  projectile/target  problem.  Since  the  mechanical  properties  in 
each  element  will  vary,  the  classically  defined  wave  velocities  are  no  longer  ap¬ 
plicable  in  the  energy  density  theory. 

Pc>-\’ndi.'nc!/  c',-;  SaT;^'acc  and  Ecdii  Fc'icci.  In  order  to  simulate  the  traction 
conditions  at  the  projectile/target  contact  surface,  it  is  necessary  to  obtain 
the  nodal  forces  in  terms  of  the  stresses  within  the  elements.  Let  these  forces 
be  denoted  as 


(28) 


They  will  act  on  the  nodes  with  lump  masses  obtained  from  the  element.  The  prin- 
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ciple  of  virtual  work  will  be  applied  to  assure  that  the  nodes  are  in  dynamic 
equilibrium.  Consider  an  arbitrary  nodal  virtual  displacement  5u®  such  that  it 
is  consistent  with  the  mechanical  constraint.  It  follows  that 


/  {oe}  {a}dV  -  /  {5u}'fT}dA  -  (6u®}{F®}  =  0 
V  '  '  A  '  '  '  ' 


in  which  {T}  is  the  traction  vector  with  normal  n  and  6e  is  the  virtual  strain 
associated  with  the  virtual  displacement  ^u.  Making  use  of  equations  (13)  and 
(22)  in  the  forms 


{iu}"^  =  {.5u®^'^[N]^ 


equation  (29)  may  be  rearranged  to  read  as 


:F®’r  +  /  [N]'  {T}dA  =  /  [B]'  f-ldV 
-  A  '  V  ' 


This  is  the  fundamental  equation  to  be  programmed  in  the  finite  element  formula¬ 
tion.  Further  simplification  of  the  volume  integral  in  equation  (32)  can  be 
made  in  view  of  the  fact  that  both  [B]  and  {j}  are  constant  over  the  triangular 
element.  For  a  unit  thickness,  it  is  obvious  that 


/  [B]'(^}dV  =  0.5[B]'{a} 


r" 


> 

;.v 

,> 

A 

'.N 

/N 


The  evaluation  of  the  surface  integral  involves  averaging  the  tractions  between 
two  nodes  and  summing  the  results. 

Unlike  the  classical  approaches  where  the  effect  of  dV/dA  is  neglected,  the 
surface  tractions  and  body  or  inertia  forces  are  no  longer  independent  in  the 
energy  density  theory.  This  effect  is  particularly  significant  at  the  interface 
where  the  projectile  comes  in  contact  with  the  target.  The  condition  to  be 
satisfied  is 


/J 


3^U  . 


T .  =  a  ■  ■  n  .  + 
T  J1  J 


'  1  i  ^  - 

3F“  MA^n’  ^ 


:34; 


with  n.  being  the  components  of  the  unit  normal  vector.  The  term  involving  the 

J 

acceleration  and  the  rate  change  of  volume  with  surface  area  is  of  the  same  order 
of  magnitude  as  the  surface  tractions  and  stresses.  The  neglect  of  (dV/dA)^  in 
equation  (34)  is  the  reason  why  the  classical  continuum  mechanics  theories  are 
not  able  to  accurately  account  for  the  transmission  of  energy  and/or  traction 
from  the  projectile  to  the  target  or  vice  versa. 

Sliding  Contact  Sux^acc.  As  the  projectile  comes  into  contact  with  the  tar¬ 
get,  sliding  of  the  adjacent  surfaces  is  permitted  and  will  be  implemented  in 
the  finite  element  formulation  following  the  technioue  presented  in  [3].  This 
necessitates  the  specification  of  the  relative  nodal  displacements  or  the  nodal 
coordinates  To  this  end,  let  the  mass  of  an  element  be  evenly  distributed 

at  the  three  nodal  points  on  which  the  nodal  forces  in  equation  (32)  are  applied. 
The  accelerations  at  these  nodes  may  be  obtained  as 


•e, 


(35) 
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where  M  stands  for  the  lumped  mass  at  the  eth  node.  The  velocity  vector  can 
thus  be  found  by  integrating  equation  (35).  Assuming  that  the  acceleration  is 
constant  for  a  small  time  increment  ^t,  then 


(36) 


*6  *  6 

Referring  to  Figure  4,  (u  }^+  and  (u  are  the  velocities  just  after  and  before 
the  current  time  t  and  At  stands  for  the  average  time  increment  about  t,  i.e., 


At 


2 


(37) 


At  incipient  impact  t  =  0,  the  initial  velocities  of  those  nodes 

in  the  projectile  which  come  in  contact  with  the  target  such  that  ^  0  for 

(At)^-.  Assuming  that  the  velocities  are  constant  over  the  time  increment  At, 
equation  (36)  can  be  integrated  to  yield  the  nodal  displacements 


t+At 


(38) 


where 


Au^;-  =  ••u^-^.(At)^+  (39) 

is,  in  fact,  the  incremental  displacement.  The  nodal  coordinates  are  then  ob¬ 
tained  from  the  displacements: 


t+At 


+  ^U 


e, 


t+At 


(40) 


-18- 


_  ’..Uf - - - Acceleration 

'  t  '  ! 


(u}  ^ - - Velocity 


‘"’t-At 


“’t+At""®'  spl acement 


Figure  4.  Time  sequence  of  displacement,  velocity  and  acceleration. 

To  reiterate,  the  nodal  velocities  are  calculated  at  the  mid-point  of  the 
time  interval  which  represents  the  average  velocities  during  the  time  increment 
At.  The  nodal  displacements  and  accelerations  are  defined  at  the  beginning  or 
the  end  of  the  time  increment.  Sliding  of  the  projectile  and  target  interface 
enables  the  readjustment  of  the  nodal  displacements,  velocities  and  accelerations 
as  a  result  of  the  interchange  of  momentum  and  kinetic  energy  between  the  projec¬ 
tile  and  target. 

Rzg^A-dding .  Because  of  the  complex  nature  of  the  projectile/target  failure 
pattern  that  changes  for  each  time  increment,  the  finite  element  grid  requires 
readjustment  as  the  system  undergoes  damage.  This  procedure  is  labor-intensive 
and  depends  on  the  experience  of  the  analyst  who  must  be  familiar  with  the  theory 
of  energy  density;  particularly,  in  grasping  the  concept  of  exchange  of  surface 
and  volume  energy.  Discussed  below  are  some  of  the  adjustments  that  must  be  con- 
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sidered  during  remeshing  of  the  finite  element  grid  pattern. 


As  the  material  between  the  projectile  and  target  is  highly  compressed, 
there  is  a  tendency  for  the  local  density  to  increase.  This  can  be  simulated  by 
judicious  redistribution  of  the  nodal  masses.  Since  the  projectile  also  suffers 
damage  in  the  form  of  fragmentation,  its  mass  will  alter  giving  rise  to  change 
of  the  projectile  momentum  and  velocity.  These  considerations  together  with 
those  outlined  in  equations  (9)  and  (10)  are  sufficient  for  analyzing  the  pro¬ 
gressive  damage  of  the  projectile/target  system. 


COMPUTER  PROGRAM 


The  computer  program  consists  of  three  main  portions.  They  are  referred  to  as 
INIXY,  PSEDA  and  PLTXY.  Table  1  shows  the  flow  chart  for  INIXY  that  deals  with 
the  initial  input  data  on  material  properties,  grid  generation  for  the  projectile/ 
target  system  and  the  relevant  parameters  that  are  required  for  carrying  out  the 
failure  analysis.  The  chart  in  Table  2  gives  an  overall  view  of  PSEDA  that  in¬ 
corporates  the  strain  energy  density  theory  for  evaluating  the  damage  and  failure 
of  all  elements.  Sliding  nodes  are  introduced  at  the  interface  where  projectile 
comes  into  contact  with  the  target  in  order  to  ensure  displacement  compatibility 
and  to  avoid  overlapping  or  interpenetration  of  material  points.  The  damage  and 
failure  pattern  for  each  time  increment  can  be  exhibited  graphically  via  the  PLTXY 
routine  shown  in  Table  3.  Contours  of  constant  dW/dV  and  dW/dA  can  also  be  plot¬ 
ted  such  that  equations  (9)  and/or  (10)  can  be  applied  to  determine  the  failure 
pattern  and  the  information  on  regridding  for  the  next  time  increment.  The  time 
increment  .‘.t  is  chosen  according  to  the  interaction  of  dynamic  disturbance  with 

the  finite  element  mesh  size.  If  '  •  denotes  the  minimum  dimension  of  the  tri- 

mi  n 

angular  element,  then  the  time  increment  should  be  chosen  to  be  smaller  than  the 


preprocess 


Node  ' -  Element 

Generation  Generat 


time  required  for  the  dynamic  disturbance  to  travel  across  5  .  ,  i.e., 

mm 

_\t  <  Min  (41) 

"o 

where  5^  is  a  normalizing  length  parameter.  The  accuracy  of  the  predictions  can 
be  improved  by  reducing  the  mesh  size  and  the  time  increment.  This  increases 
the  labor  of  regridding  and  computer  time  in  addition  to  the  readjustments  that 
must  be  made  on  mass  redistribution  and  momentum  change.  Refer  to  Tables  4,  5 
and  6  for  a  more  detailed  account  of  the  computing  procedures  involving  INIXY, 
PSEDA  and  PLTXY. 

IMPACT  OF  PROJECTILE  ON  TARGET  PLATE 

The  computer  program  described  previously  will  be  referred  to  as  the  Plane 
Strain  Energy  Density  Analysis  (PSEDA)  code  and  it  will  be  applied  to  solve  the 
problem  of  a  projectile  impacting  a  target  plate.  Analyzed  will  be  the  sequence 
of  damage  pattern  as  the  projectile  penetrates  through  the  target  in  an  incre¬ 
mental  manner  without  making  any  assumptions  on  the  failure  modes  and  the  local 

strain  rates.  The  material  data  bank  has  the  capacity  of  covering  strain  rates 
-4  6  -1 

from  10  to  10  sec  which  is  more  than  adequate  for  the  example  to  be  con¬ 
sidered. 


Sijitcir^  GL’ometUf.  A  tungsten  projectile  with  dimensions  Dp  =  1.4  cm  and  Lp 
4.7  cm  is  traveling  at  a  velocity  of  1,200  m/sec  and  comes  into  normal  contact 


with  a  plate  with  dimensions  h^  =  2  cm  and  D^  =  10  cm  as  illustrated  in  Figure  5 
Because  of  symmetry,  only  one-half  of  the  problem  needs  to  be  discretized  by  ap¬ 
plication  of  triangular  elements.  The  initial  grid  pattern  is  displayed  in  Fig¬ 
ure  6  in  which  64  elements  and  52  nodes  are  employed.  Smaller  elements  are  de- 


»  •>  ->->  ■> 


Table  6.  Block  diagram  for  computing  algorithm  of  PLTXY. 
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picted  near  the  contacting  end  with  a  reentrant-corner  where  the  stress  and  en¬ 
ergy  are  concentrated  and  failure  is  most  likely  to  initiate.  The  locations  of 
each  element  and  node  are  referenced  to  the  x-  and  y-coordinate.  An  enlarged 
view  of  the  region  near  the  contact  end  is  given  in  Figure  7  which  also  provides 
information  on  the  numbering  of  the  nodes  and  elements.  As  damage  proceeds,  the 
grid  pattern  in  Figure  6  will  be  modified  accordingly. 

-3  -4  -1 

Bcxic  'lafottaf.  Under  normal  strain  rates,  say  10  or  10  sec  ,  the  stress 
and  strain  response  of  most  metals  will  first  rise  sharply  and  then  level  off  un¬ 
til  fracture  occurs  at  f  in  Figure  8  at  which  point  the  uniaxial  test  terminates. 
Such  a  behavior,  however,  will  not  be  valid  for  all  points  in  the  projectile/tar¬ 
get  system  after  the  initial  impact  because  the  local  strain  rates  will  not  only 
be  higher  but  they  will  vary  from  one  location  to  another.  Therefore,  the  base 
material  property  can  be  used  only  during  the  early  onset  of  impact. 

As  a  reference,  Table  7  gives  the  conventional  mechanical  properties  of  the 

tungsten  projectile  and  the  4340  steel  target  plate.  The  symbols  E>  c  and  zr 

ys  1 

stand,  respectively,  for  the  Young's  modulus,  yield  strain  and  final  strain.  The 
coefficients  =  and  y  correspond  to  those  used  in  the  uniaxial  stress  and  strain 
relation 

7  -  1]^  (42) 
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The  resulting  curves  are  given  in  Figure  9.  Note  that  the  tungsten  has  a  higher 

yield  strength  a  and  fracture  toughness  related  to  (dW/dV)  or  area  under  the 
ys  c 

true  stress  and  true  strain  curve  than  the  4340  steel.  This  is  further  illus¬ 
trated  in  Figure  10  where  *  tends  to  drop  as  (dW/dV)  increases  for  both  the 

ys  c 

projectile  and  target  material.  Such  a  trade-off  is  well-known  for  metals  and 


-30- 


Stress  .)  X  10  (N/in' 


Figure  10.  Relation  between  yield  strength  and  critical  strain 
energy  density  function  for  projecti 1  e/target  system 
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Figure  11.  Initial  geometry  of  projectile/target  system 
at  t  =  0.0  usec. 


There  is  no  formal  procedure  for  optimizing  the  finite  element  mesh  size  with 
the  corresponding  time  step.  Any  variance  will  alter  the  energy  dissipation 
history  or  the  time  dependency  of  the  H-function  [11,12]  that  determines  the 
degree  of  irreversibility  or  dissipative  character  of  the  projectile/target  sys¬ 
tem.  A  thorough  understanding  of  the  combined  influence  of  mesh  size  and  time 
step  on  the  damage  of  a  physical  system  in  general,  however,  is  a  topic  beyond 
the  scope  of  this  investigation. 

The  contours  of  dW/dA  on  the  damage  plane  and  dW/dV  are  first  obtained  for 

t  =  0.075  usec.  They  are  displayed  in  Figures  12  and  13.  To  be  kept  in  mind 

is  that  the  critical  values  (dW/dA)^  and  (dW/dV)^  will  change  from  location  to 

location  because  the  stress  and  strain  response  in  each  element  is  different. 

For  instance,  Figure  12  shows  that  the  intensity  dW/dA  for  contours  no.  1  to 

3  3 

no.  4  increases  from  0.86  x  10  N/m  to  27.95  x  10  N/m.  It  is  not  obvious  where 
failure  will  first  initiate  unless  the  (dW/dA)^  values  at  these  locations  are 
carefully  examined.  According  to  the  condition  in  equation  (10),  failure  by 
localized  spall  occurred  first  at  the  front  surface  of  the  target  directly  under¬ 
neath  the  projectile  corner.  This  justifies  the  removal  of  a  semi-circular  block 
of  material  whose  boundary  coincides  with  the  surface  along  which  (dW/dV)^ 

3 

=  14.88  X  10  N/m.  Even  though  the  dW/dA  values  of  contours  no.  3  and  no.  4  are 
higher,  no  failure  is  predicted  at  these  locations  because  the  corresponding 
(dW/dA)^  values  are  higher  than  19.21  x  10^  N/m  and  27.95  x  10^  N/m.  Figure  13 
shows  that  the  magnitude  of  the  dW/dV  contours  remained  below  the  critical  value. 
A  corner  cavity  is  thus  modelled  in  the  target  after  a  time  lapse  of  0.075  usec. 


It  is  analogous  to  the  entropy  function  in  classical  thermodynamics  when  energy 
is  dissipated  all  in  the  form  of  heat.  In  general,  H  acquires  a  much  broader 
interpretation  [11,12]  of  order  and  disorder  as  disturbed  by  loading  and  can  be 
either  negative  or  positive. 
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Contours  of  constant  dW/dA  at  t  =  0.075  ..sec. 
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Figure  13.  Contours  of  constant  dW/dV  at  t  =  0.075  usec 


Because  the  cavity  is  sufficiently  small,  the  same  number  of  64  elements  and  52 
nodes  can  be  used  for  the  next  time  increment  except  that  the  elements  next  to 
the  cavity  have  been  distorted  as  indicated  in  Figure  15.  Calculations  made  at 
t  =  0.120  usec  yield  another  set  of  dW/dA  and  dW/dV  contours  and  they  revealed 
increasing  dynamic  disturbances  in  both  the  projectile  and  target  as  shown  in  Fig 
ures  16  and  17.  Again,  all  dW/dV  contours  are  below  critical  while  the  critical 

3 

(dW/dV)^  =  16.04  X  10  N/m  contour  prevailed  in  a  region  circling  the  cavity  in 
Figure  14  or  15.  An  increase  in  the  cavity  size  is  thus  predicted  as  shown  in 
Figure  18.  Remeshing  of  the  finite  element  grid  pattern  is  carried  out  in  Fig¬ 
ure  19  such  that  the  elements  around  the  cavity  periphery  are  made  smaller  for 
better  accuracy.  This  results  in  an  increase  of  element  number  from  64  to  67 
and  nodal  points  from  52  to  56. 

As  the  material  in  the  target  under  the  projectile  corner  is  being  removed, 
energy  transmission  begins  to  shift  towards  the  center  of  the  target.  For  t 
=  0.175  usec,  the  values  of  (dW/dV)^  near  the  target  surface  under  and  away  from 

3  3 

the  projectile  are,  respectively,  15.53  x  10  N/m  and  14.88  x  10  N/m.  Accord¬ 
ing  to  the  dW/dA  contours  displayed  in  Figure  20,  further  fragmentation  of  the 
target  surface  is  predicted.  The  corresponding  dW/dV  contours  no.  1,  2  and  3  in 
the  same  neighborhood  are  of  the  order  of  10^  N/m^  which  is  still  one  order  of 
magnitude  lower  than  the  critical  value.  Surface  energy  density  remains  as  the 
dominant  cause  of  failure.  Figure  22  shows  the  fragmented  shape  of  the  target 
surface  and  the  corresponding  grid  pattern  is  given  in  Figure  23.  No  surface 
separation  has  yet  occurred  within  the  projectile.  Damage  continues  to  occur  as 
t  advances  from  0.175  usec  to  0.240  usec.  The  change  in  the  dW/dA  contours  is 
small.  This  can  be  seen  by  comparing  the  results  in  Figure  24  with  those  in  Fig¬ 
ure  20.  Higher  values  of  dW/dV  around  the  fractured  region  of  the  target  did 
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occur  in  a  small  element  as  indicated  in  Figure  25  with  (dW/dV)^  =  7.29  x  10  M/m 
A  slight  alteration  of  the  fragmented  zone  is  made  in  Figure  26  although  its  dif¬ 
ference  from  Figure  22  is  not  noticeable  at  the  scale  level  presented.  Refer  to 
Figure  27  for  the  remeshed  grid  pattern  after  0.240  usec. 


Since  the  dynamic  disturbances  propagate  through  the  system  at  such  a  rapid 
rate,  the  geometry  of  the  projectile/target  system  remained  virtually  unchanged. 
Exhibited  in  Figure  28  are  the  dW/dA  contours  in  both  the  projectile  and  target 
for  t  =  0.345  usec.  At  the  projectile  corner,  dW/dA  attained  the  values  of  41.52 
X  10  N/m  and  60.56  x  10“^  N/m  which  have  exceeded  the  critical  surface  energy 

3 

density  of  38.99  x  10  N/m.  Fracture  of  the  corner  element  is  thus  assumed  to 
take  place.  Additional  damage  also  prevails  near  the  already  fragmented  surface 
of  the  target  where  (dW/dA)^  has  reached  the  values  of  10.9  x  10^  N/m  and  16.2 
X  10  N/m  as  shown  in  Figure  28.  Similarly,  the  dW/dV  contours  in  Figure  29  show 
that  (dW/dV)  =  179.6  x  10^  M/m-  at  the  projectile  corner  and  (dlJ/dV)  =  80.34 

C 

X  10  N/m^  next  to  the  cavity  boundary  in  the  target  have  also  been  exceeded. 

The  resulting  failure  pattern  is  summarized  schematically  in  Figure  30  while  the 
corresponding  grid  pattern  is  presented  in  Figure  31.  From  here  on,  damage  in 
the  projectile  and  target  tends  to  spread  more  widely  as  illustrated  by  the  en¬ 
ergy  density  contours  in  Figures  32  and  33.  The  computer  output  reveals  that  at 
t  =  0.480  usec,  the  (dW/dV)^  =  11.6  x  10  N/m  critical  boundary  extends  between 
contours  no.  1  and  no.  2  in  Figure  32  from  the  target  surface  to  its  mid-portion 
dislocating  and  fracturing  a  large  portion  of  the  material.  Figure  33  shows  that 
dW/dV  in  several  of  the  protruded  regions  has  also  been  severely  intensified  be¬ 
yond  the  critical  value  of  (dW/dV)^  =  85.98  x  10^  N/m-.  The  material  in  these 
regions  is  fragmented  into  small  pieces  and  can  be  removed  by  redistribution  of 
the  mass  onto  the  neighboring  nodal  points.  A  schematic  of  the  project! le/taraet 
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Figure  28.  Contours  of  constant  dW/dA  at  t  =  0.345  usec. 
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Figure  29.  Contours  of  constant  dW/dV  at  t  =  0.345  usec 

-56- 


N  ■>  ->  --  ^ 


1.6  0.0 


H - 1 — 

1.6  3.2 

X  (cm) 


tTs 


Figure  30.  Project! 1 e/ target  damage  pattern  after  0.345  .sec 
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Figure  32.  Contours  of  constant  dW/dA  at  t  =  0.480  usec 


failure  pattern  after  0.480  ysec  is  shown  in  Figure  34.  Three  separate  frag¬ 
ments  are  indicated  which  are  in  loose  contact  with  the  target.  A  gap  follow- 
ing  the  (dW/dV)^  =  11.6  x  10  N/m  critical  contour  is  modelled  in  Figure  35 
that  contains  114  elements  and  97  nodes.  The  largest  block  of  loose  material  is 
near  the  center  and  is  barely  in  contact  with  the  projectile.  The  damage  pattern 
at  t  =  0.501  usec  of  the  projectile/target  system  is  virtually  unchanged  from 
that  at  t  =  0.480  usec.  What  has  altered  slightly  is  only  the  partial  disinte¬ 
gration  of  one  of  the  fragments  as  illustrated  in  both  Figures  36  and  37.  This 
minor  correction  is  made  in  Figures  38  while  Figure  39  gives  the  new  grid  pat¬ 
tern  for  the  next  time  increment  of  computation.  The  results  for  t  =  0.562  usec 
are  displayed  in  Figures  40  and  41.  As  the  dynamic  disturbances  have  now  been 
well  propagated  into  the  system,  the  surface  energy  density  dW/dA  begins  to  be 
less  dominant  and  remains  below  the  damage  threshold  shown  in  Figure  40  while  the 
volume  energy  density  dW/dV  starts  to  increase  more  readily.  Further  disintegra¬ 
tion  of  the  fragments  occurred  as  the  volume  energy  density  intensified  substan¬ 
tially  beyond  (dW/dV)^  =  85.98  x  10®  N/m^,  Figure  41.  Modified  in  Figures  42 
and  43  are  the  fragmented  configuration  after  0.562  usec. 

Damage  continues  to  accumulate  as  t  is  increased  to  0.663  usec  due  to  the  in¬ 
crease  in  the  local  values  of  dW/dA  and  dW/dV.  The  sites  of  failure  owing  to  ex¬ 
cessive  surface  energy  are  identified  in  Figure  44  with  the  different  values  of 
(dW/dA)^  in  the  target  and  projectile.  Figure  45  gives  the  corresponding  criti¬ 
cal  values  of  dW/dV.  The  large  fragment  that  was  originally  in  contact  with  the 
projectile  has  now  been  broken  into  two  resulting  in  a  total  of  four  pieces  and 
leaving  a  space  between  the  projectile  and  target.  This  is  shown  in  Figure  46 
with  the  corresponding  grid  pattern  in  Figure  47  that  employed' 100  elements  and 
89  nodes.  Multiple  fragmentation  after  0.633  usec  created  a  wide  space  such 
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Figure  37.  Contours  of  constant  dW/dV  at  t  =  0.501  ..sec 
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Figure  39.  Grid  pattern  of  projectile/target  system  after  0.501  usec. 
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Figure  44.  Contours  of  constant  dW/dA  at  t  =  0.663  usec. 
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that  the  projectile  is  not  in  contact  with  the  target.  The  local  surface  and 
volume  energy  density  as  shown  in  Figures  48  and  49  continue  to  elevate  leading 
to  an  increase  in  the  cavity  size  and  further  disintegration  of  the  fragments. 

The  smaller  debris  has  little  or  no  influence  on  the  dynamic  behavior  of  the 
system  and  can  be  removed.  This  leads  to  the  configuration  in  Figure  50  whose 
grid  pattern  is  given  in  Figure  51.  A  glance  at  the  results  in  Figures  52  to  55 
show  that  almost  the  entire  projectile  is  now  under  stress  or  strain.  However, 
change  in  the  failure  pattern  during  the  time  interval  of  t  =  0.764  usec  and 
0.869  usec  involves  only  the  disintegration  of  one  fragment  near  the  center. 

There  remains  only  two  loose  fragments  as  indicated  in  Figures  54  and  55.  Even 
though  the  projectile  is  not  in  direct  contact  with  the  target  at  t  =  0.924  usec, 
damage  continues  to  occur  because  of  the  increase  in  dW/dA  and  dW/dV.  Separation 
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of  a  layer  of  material  or  fracture  is  predicted  near  the  projectile  end  in  Figures 
56  and  57  because  both  the  surface  and  volume  energy  density  have  exceeded  the 
respective  critical  values  of  25.86  x  10  N/m  and  180.5  x  10  N/m^.  The  volume 
energy  density  in  two  fragments  is  also  much  greater  than  (dW/dV)^  =  48  x  10^ 

N/m^.  At  this  instance,  all  fragments  have  disintegrated  leaving  a  large  cavity 
between  the  projectile  and  target.  This  configuration  is  modelled  in  Figures  58 
and  59. 
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Since  the  projectile  is  still  not  in  contact  with  the  target,  there  is  time 
for  the  surface  energy  density  near  the  cavity  boundary  to  decrease  below  the 
critical  value.  Figure  60.  The  volume  energy  density  inside  the  surface  given 
in  Figure  61,  however,  remains  large  and  causes  further  fracture.  A  zig-zag 
fracture  surface  is  predicted  at  t  =  1.024  usec.  Refer  to  Figure  62  and  the  grid 
pattern  in  Figure  63.  Again,  the  surface  energy  density  near  the  fractured  area 
in  the  target  remained  below  critical  while  only  that  in  the  projectile  as  shown 
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Figure  59.  Projectile  and  target  grid  pattern  at  t  =  0.924  usec 
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in  Figure  64  exceeded  (dW/dV)^.  The  volume  energy  density  contours  at  t  =  1.224 
ysec  in  Figure  65,  however,  have  increased  beyond  (dW/dV)^  at  several  locations. 
This  results  in  additional  fragmentation  as  displayed  in  Figures  66  and  67.  Fur¬ 
ther  calculations  made  for  t  =  1.245  usec  revealed  no  significant  change  except 
for  the  disintegration  of  the  fragment  that  is  in  point  contact  with  the  frac¬ 
tured  surface. 

The  cavity  depth  in  the  target  at  t  =  1.259  usec  is  approximately  one-half 
of  the  total  thickness.  Even  though  dW/dA  is  relatively  low  as  indicated  in  Fig¬ 
ure  68,  the  volume  energy  density  dW/dV  in  several  places  is  still  above  critical. 
Failure  by  fracture  occurs  at  (dW/dV)^  =  32.22  x  10^  N/m^  in  the  target  and 
(dW/dV)^  =  288.37  x  10^  N/m^  in  the  projectile  as  shown  in  Figure  69.  A  large 
portion  of  the  material  breaks  away  from  the  projectile  which  continues  to  travel 
toward  the  target  as  a  rigid  body.  This  is  illustrated  in  Figure  70.  Figure  71 
displays  the  revised  grid  pattern  for  the  next  time  increment  of  calculation. 

As  the  projectile  comes  in  contact  with  the  trapped  fragment  at  the  interface, 
dW/dA  remained  relatively  low  everywhere.  This  is  illustrated  in  Figure  72  for 
t  =  1.264  usec.  The  fragment  is  crushed  as  a  result  of  the  high  volume  energy 
density  surpassing  the  critical  value  of  129.46  x  10^  N/m  indicated  in  Figure  73. 

A  gap  between  the  projectile  and  target  is  again  created  and  shown  in  Figure  74. 
The  corresponding  finite  element  grid  pattern  in  Figure  75  is  now  simplified  and 
contains  only  61  elements  and  49  nodes.  Once  the  projectile  catches  up  to  the 
cavity  in  Figure  76,  the  local  surface  energy  density  intensifies  and  rises  above 
(dW/dV)^  =  29.59  x  10  N/m.  The  local  volume  energy  density  contours  in  Figure 
77  also  attain  high  values  but  they  are  localized.  It  can  be  easily  seen  from 
the  computer  output  that  there  exists  a  (dW/dA)^  critical  path  extending  from 
the  reentrant  corner  of  the  cavity  in  Figures  76  or  77  to  the  back  side  of  the 

target  plate.  This  event  is  normally  referred  to  as  "plugging"  where  a  block  of 
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Figure  74.  Damage  pattern  at  t  =  1.264  usec 
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tarcjet  material  is  dislodged  and  leaves  with  the  projectile.  The  final  damage 
configuration  of  the  projectile/target  system  is  shown  schematically  in  Figure 
78. 


CONCLUDING  REMARKS 

Prediction  on  the  progressive  failure  of  a  projectile/target  system  ading 
up  to  plug  formation  has  been  presented  by  application  of  the  energy  density 
theory.  No  a  priori  assumption  were  made  on  the  constitutive  relations  of  the 
materials  that  are  not  only  highly  stressed  but  their  strain  rates  change  with 
location  and  time.  Each  time  interval  was  carefully  chosen  such  that  a  complete 
sequence  of  the  damage  pattern  is  obtained  capturing  all  of  the  failure  modes 
including  permanent  deformation,  fragmentation,  fracture  and  plug  formation. 

This  corresponds  specifically  to  a  tungsten  projectile  with  an  aspect  ratio  of 
1:3.36  impacting  on  a  4340  steel  target  with  an  aspect  ratio  of  1:5  at  an  ini¬ 
tial  velocity  of  1,200  m/sec.  The  sequence  of  failure  modes  depend  sensitively 
on  the  combination  of  material  type,  projectile/ target  geometry  and  size  and  the 
initial  impact  velocity.  For  the  problem  at  hand,  the  following  main  sequence  of 
events  may  be  summarized: 

(1)  Time  t  up  to  0.480  usec:  Owing  to  the  strain  and/or  energy  concentra¬ 
tion  around  the  projectile  periphery,  a  cavity  is  created  in  the  tar¬ 
get.  This  shifts  the  load  transfer  to  the  target  material  at  the  center 
which  is  still  in  contact  with  the  projectile. 

(2)  Time  from  t  =  0.480  to  0.663  usec:  Because  of  the  high  load  transfer 


over  a  small  portion  of  the  contacting  target  material,  fragmentation 
and  disintegration  of  debris  occurred  leaving  a  large  void  or  cloud  of 
small  fragments  between  the  projectile  and  target. 


(3)  Time  from  t  =  0.663  to  0.924  ysec:  The  projectile  although  strained 


and  disturbed  dynamically  travels  essentially  as  a  rigid  body  to  close 
the  gap  created  by  fragmentation. 

(4)  Time  from  t  =  0.924  to  1.224  psec:  As  the  projectile  catches  up  to  the 
partially  damaged  target  that  is  still  highly  strained,  the  impact 
again  results  in  local  fragmentation  which  is  not  as  severe  as  that 
which  occurred  during  the  time  interval  from  t  =  0.480  to  0.663  usec. 

(5)  Time  from  t  =  1.224  to  1.265  usec:  During  this  time  interval,  all 
fragmentation  is  completed  and  the  remaining  energy  from  the  projectile 
is  concentrated  along  a  path  that  outlines  a  plug.  Its  formation  repre¬ 
sents  the  onset  of  complete  penetration. 

P^.nc  traticn  Tunc  and  Ex.Lt  Vciocitij.  An  estimate  on  the  time  for  the  projec¬ 
tile  to  penetrate  completely  through  the  target  can  be  made  by  considering  both 
the  time  consumed  in  the  damage  progress  in  addition  to  that  required  for  the 
projectile  to  travel  as  a  rigid  body.  The  result  must  necessarily  be  greater 
than  16.67  usec  which  corresponds  to  the  time  of  the  original  projectile  travel¬ 
ling  at  1,200  m/sec  through  a  distance  of  2  cm,  the  thickness  of  the  target  plate 
Table  8  summarizes  the  time  increments  elapsed -for  a  total  of  18  steps  such  that 
the  rigid  body  time  for  the  projectile  is  also  included  beyond  step  no.  19.  Note 
that  the  total  time  up  to  the  onset  of  plug  formation  is  20.085  asec.  The  addi¬ 
tional  time  required  for  the  projectile  and  plug  to  leave  the  target  is  11.249 
usec.  This  gives  a  value  of  31.325  usec  at  which  time  the  plug  exits  with  a  ve¬ 
locity  of  802.3  m/sec.  The  time  history  of  the  projectile  velocity  is  displayed 
graphically  in  Figure  79. 
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Table  8.  Time  increments  leading  up  to  plugging  in  usec. 


ep  No. 

T ime  Increment 

Rigid  Body 

Total  Time 

1 

0.075 

0.075 

2 

0.045  (0.120)* 

0.120 

3 

0.055  (0.175) 

0.175 

4 

0.065  (0.240) 

0.240 

5 

0.105  (0.345) 

0.345 

6 

0.135  (0.480) 

0.480 

7 

0.021  (0.501) 

0.501 

8 

0.061  (0.562) 

0.562 

9 

0.101  (0.663) 

0.663 

10 

0.101  (0.764) 

0.003 

0.767 

11 

0.105  (0.869) 

0.008 

0.880 

12 

0.055  (0.924) 

0.006 

0.941 

13 

0.100  (1.024) 

0.014 

1.056 

14 

0.200  (1.224) 

0.348 

1.604 

15 

0.021  (1.245) 

0.764 

2.389 

16 

0.014  (1.259) 

3.619 

6.023 

17 

0.005  (1.264)** 

9.644 

15.671 

18 

0.001  (1.265) 

4.912 

20.085 

Time  increment  excluding  rigid  body  motion. 

Time  step  not  presented  in  text. _ 

C!:an{(C  n:  U'cai  Strain  Rate.  As  mentioned  earlier,  each  element  in  the  pro¬ 
jectile/target  system  experiences  a  different  strain  rate.  The  local  values  de¬ 
pend  on  the  area  over  which  the  averages  are  taken.  As  an  example,  consider  the 
grid  pattern  in  Figure  80  for  t  =  1.265  usec  that  corresponds  to  the  instant  of 
plug  initiation.  The  values  of  n,  and  on  the  damage  plane  are  averaged  over 
the  elements  no.  220,  221, — ,  247  covering  the  region  where  failure  is  antici¬ 
pated.  Their  values  are  given  in  Table  9  and  plotted  in  Figure  81  for  some  typi- 


Projectile  velocity  x  10  (m/sec) 


cal  elements.  Refer  to  Table  8  for  the  definition  of  time  elapsed.  The  dif¬ 
ference  in  the  dynamic  stress  and  strain  response  in  the  elements  near  the  pro¬ 
spective  site  of  plug  formation  is  obvious.  The  highest  strain  rate  of  10^  sec' 
occurred  in  element  fio.  220  and  227  whose  sizes  are  large  in  comparison  with  the 
local  region  next  to  the  so-called  adiabatic  shear  band  which  has  been  observed 
experimentally  [8,9]. 

A  much  higher  local  strain  rate  up  to  10^  sec'"'  would  result  if  the  results 
are  averaged  over  elements  with  linear  dimensions  of  the  same  order  as  the  shear 
band.  Inaccuracy  is  also  anticipated  in  the  region  near  the  backside  of  the  tar 
get  where  local  spallation  may  occur.  This  can  be  determined  by  resulting  to 
smaller  elements,  a  task  that  will  be  left  for  future  investigation. 

Ti:o ‘uriac 'Mcccatitcac  A  salient  feature  of  the  energy  density  theory 

is  that  thenmal/mechanical  effects  are  automatically  considered  as  inherent  in 
the  physical  damage  process.  Because  unloading  is  calculated  rather  than  preas¬ 
sumed,  dissipation  can  be  determined  for  each  time  increment  or  as  a  function  of 
time.  Without  loss  in  generality,  the  volume  energy  density  dW/dV  in  equation 
(1)  or  (26)  can  be  divided  into  two  parts  denoted  as  (dW/dV)p  and  (dW/dV)  ,  i.e. 


dV  '  MV^p 


in  which  P  =  (dW/dV)p  is  the  dissipation  energy  density  and  A  =  (dU/dV)  the 
available  energy  density.  Their  sum  represents  the  area  under  p  as  referred  to 
the  true  stress  and  true  strain  diagram  in  Figure  82.  Once  P  is  known,  the  tem¬ 
perature  in  each  element  can  be  found  from  [6,7] 


True  strain  c. 

Figure  82.  Dissipated  and  available  energy  density  diagram. 


where  and  are  the  equivalent  uniaxial  stress  and  strain  on  the  damage  plane, 

4 

The  quantity  AP/Ar.,  in  equation  (44)  can  be  associated  with  the  latent  heat  in 
classical  thermodynamics  for  evaluating  change  in  material  microstructure  or 
phase  transformation,  say  the  melting  of  solid  or  the  evaporation  of  liquid.  If 
dissipation  takes  place  all  in  the  form  of  heat,  then  P  0  and  '  T.  It  fol¬ 
lows  that 


^  _  T  ^ 
.tV  -  tV 


Here,  tQ  is  the  heat  required  per  unit  mass  of  substance  in  changing  the  thermo¬ 
dynamic  state  from  T  to  T+;.T  and  .tV  is  the  change  of  specific  volume.  The  corre¬ 
sponding  change  in  entropy  is  ,',S.  Equation  (45)  can  be  used  to  calculate  the 


I 
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latent  heat  at  a  given  rate  and  temperature  in  the  4340  steel  near  the  site 
where  the  adiabatic  shear  band  was  formed  when  plugging  occurred.  A  check  can 
thus  be  available  to  verify  whether  the  martensite  microstructure  indeed  changes 
to  the  austenitic  microstt^ucture  which  was  observed  experimentally.  Such  a 
transformation  is  dc^Kndcnt ,  an  effect  that  cannot  be  accounted  for  in 

classical  thermodynamics. 

There  is  no  difficulty  to  include  phase  transformation  in  the  energy  density 
theory.  Projectile/target  failure  behavior  in  the  hypervelocity  range,  say 
9,000  m/sec,  can  also  be  analyzed  where  the  target  material  under  the  projectile 
can  transform  into  liquid  and/or  gas.  The  damaged  zone  of  the  target  may  change 
drastically  from  the  original  configuration.  Preliminary  work  on  hypervelocity 
impact  has  been  initiated  at  the  Institute  of  Fracture  and  Solid  Mechanics  at 
Lehigh  University. 
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